#include<iostream>
#include<cstring>
#include<vector>
using namespace std;
const int N=5e2+10;
int n,m,s,d;
int g[N][N],dis[N],vis[N];
int p[N];
long long cnt=0,sum=0;
vector<int>res,tem;
void dij()
{
	memset(dis,0x3f,sizeof dis);
//	for(int i=0;i<n;i++)
//	dis[i]=g[s][i];
	dis[s]=0;
	for(int i=0;i<n;i++)
	{
		int cur=-1;
		for(int j=0;j<n;j++)
		{
			if(!vis[j]&&cur==-1||dis[cur]>dis[j])
			cur=j;
		}
		vis[cur]=1;
		for(int j=0;j<n;j++)
		{
			dis[j]=min(dis[cur]+g[cur][j],dis[j]);
		}
	}
}
void dfs(int u,long long v,long long s)
{
	if(v>dis[d])
	return ;
	if(u==d&&v==dis[d])
	{
		cnt++;
		if(s>sum)
		sum=s,res=tem;
		return ;
	}
	for(int i=0;i<n;i++)
	{
		if(!vis[i])
		{
			vis[i]=1;
			tem.push_back(i);
			dfs(i,v+g[u][i],s+p[i]);
			tem.pop_back();
			vis[i]=0;
		}
	}
}
int main()
{
	memset(g,0x3f,sizeof g);
	cin>>n>>m>>s>>d;
	for(int i=0;i<n;i++)
	cin>>p[i];
	for(int i=0;i<m;i++)
	{
		int u,v,w;
		cin>>u>>v>>w;
	    g[v][u]=g[u][v]=min(g[u][v],w);
	}
	dij();
	memset(vis,0,sizeof vis);
	dfs(s,0,p[s]);
	cout<<cnt<<' '<<sum<<endl;
	cout<<s;
	for(int i=0;i<res.size();i++)
	cout<<' '<<res[i];
	return 0;
}
